Electronic structure and the glass transition in pnictide and chalcogenide 
semiconductor alloys. Part I: The formation of the ppa- network. 
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Semiconductor glasses exhibit many unique optical and electronic anomalies. We have put forth 
a semi-phenomenological scenario (J. Chem. Phys. 132, 044508 (2010)) in which several of these 
anomalies arise from deep midgap electronic states residing on high-strain regions intrinsic to the 
activated transport above the glass transition. Here we demonstrate at the molecular level how this 
scenario is realized in an important class of semiconductor glasses, namely chalcogen and pnictogen 
containing alloys. Both the glass itself and the intrinsic electronic midgap states emerge as a result of 
the formation of a network composed of u-bonded atomic p-orbitals that are only weakly hybridized. 
Despite a large number of weak bonds, these ppa- networks are stable with respect to competing types 
of bonding, while exhibiting a high degree of structural degeneracy. The stability is rationalized 
with the help of a hereby proposed structural model, by which ppcr-networks are symmetry-broken 
and distorted versions of a high symmetry structure. The latter structure exhibits exact octahedral 
coordination and is fully covalently-bonded. The present approach provides a microscopic route to 
a fully consistent description of the electronic and structural excitations in vitreous semiconductors. 



I. INTRODUCTION 



The electronic and structural excitations in amor- 
phous semiconductors, and the interplay of these exci- 
tations, have evaded a self-consistent first-principles de- 
scription for decades. Amorphous semiconductors^— are 
important both in applications, e.g., as phase change 
materials?^ and from the basic viewpoint. The elec- 
tronic structure in a disordered lattice is fundamentally 
different from the venerable Bloch picture of continuous 
bands of allowed states separated by strictly forbidden 
gaps, as would be applicable in periodic solids. Although 
a result of multiple electron scattering, the presence of 
such forbidden gaps in periodic solids is also consistent 
with their energetic stability in that such gaps usually 
imply stabilized occupied orbitals. In contrast, in a dis- 
ordered lattice, strict gaps in the density of states are 
not a priori allowed. Yet the electronic orbitals can 
still be subdivided in two relatively distinct classes^— 
(a) extended states, in which electrons move as well- 
defined wave-packets, and (b) localized states, whose 
density decays nearly exponentially away from mobil- 
ity bands. Despite the many successes in applying these 
ideas semi-phenomenologically, developing a first princi- 
ples description, in which a realistically stable aperiodic 
lattice and mobility gaps emerge self-consist entlyj 1 1 ' 12 
has been difficult^ Further compounding this difficulty, 
several electronic and optical peculiarities of disordered 
semiconductors indicate that there are effects of disor- 
der beyond those generically expected of a mechanically 
stable disordered lattice. 

We have argued 14 that, instead, it is not the aperiodic- 
ity alone, but the intrinsic metastability of semiconductor 
glasses that leads to many unique properties of these dis- 
ordered solids. A glass is made by thermally quenching a 
supercooled liquid at a rate exceeding the typical liquid 
relaxation rate. A supercooled liquid can be thought of 



as an aperiodic crystal characterized by myriad low free 
energy configurations that are nearly degenerate. Molec- 
ular motions in the liquid occur via local activated tran- 
sitions between these low free energy configurations? 1 ^— 
which are accompanied by the creation of high-strain in- 
terfacial regions separating the configurations. The in- 
terfaces are intrinsic to the activated dynamics; their 
concentration, for a given quenched glass, depends (loga- 
rithmically slow) only on the time scale of the glass tran- 
sition, i.e., the quench speed. According to the Ran- 
dom First Order Transition (RFOT) theory of the glass 
transition,— ~— this concentration just above the glass 
transition temperature on one hour scale, T g 



cally about £" 



10 cm , within an order of magni- 



tude or so, depending on the specific substance. The pa- 
rameter £ denotes the cooperativity length for activated 
reconfigurations . 



Our results indicate^ that the high-strain regions 
that form when the amorphous material is made may 
host deep midgap electronic states of topological ori- 
gin, which are centered on over- and under-coordinated 
atoms. These states share several characteristics with the 
midgap electronic states in trans-polyacetylene, which 
are centered on defects in the perfect alternation pat- 
tern of the double and single bonds along the polymer 
chains. 19 Using a semi-phenomenological, coarse-grained 
Hamiltonian, we have established the spatial and charge 
characteristics of the interface-based midgap states in 
non-polymeric glasses. We further concluded, based on 
the internal consistency of the description, that the states 
should be present only in a limited class of glasses that 
satisfy the following requirements: The bonding should 
exhibit inhomogeneous saturation so that the transfer in- 
tegrals t in the electronic effective tight-binding Hamil- 
tonian H = V, (,<■]<■ , + J2{ij) tij c l c j snou ld uniformly 
exhibit spatial variation. Nevertheless, the magnitude of 
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the variation should be modest: 14 

1*7*1 £0.5, (1) 

where t and t' denote the upper and lower limits of the 
variation range. More detailed estimates 2 ^ indicate that 
the lower limit on the t' jt ratio is probably smaller, i.e., 
0.3 or so. Finally, the spatial variation Se in electroneg- 
ativity should not be too large 

\5e\<\t-t'\, (2) 

thus implying the material is a semiconductor, since 
the transfer integral t is at most a few eV. Of 
non-polymeric materials, only certain chalcogen- and 
pnictogen-containing glasses appear to satisfy all of 
these requirements. On the other hand, these amor- 
phous arsenic chalcogenides and similar materials do 
indeed exhibit several electronic and optical anoma- 
lies that could be accounted for by the interface-based 
states, in a unified fashion^ These anomalies include 
light-induced electron spin resonance (ESR) and midgap 
absorption ) 21 i 22 two types of photoluminescence^ and 
field-induced ultrasonic attenuation^ Thus general ar- 
guments, on the one hand, and observation, on the other, 
seem to converge on the uniqueness of chalcogenide and 
pnictide glasses with regard to their potential ability to 
host topological midgap states. Despite this remark- 
able convergence, the currently available evidence for the 
unique interplay between electronic excitations and the 
metastability in those glasses must be regarded as cir- 
cumstantial. 

The purpose of the present effort is to test the conclu- 
sions of the semi-phenomenological analysis from Rcf. 14 
directly, based on the local chemistry specific to amor- 
phous chalcogenide and pnictide alloys. Our basic hy- 
pothesis for the origin of the topological midgap states 
in the semiconductor glasses is that these glasses repre- 
sent aperiodic networks of cr-bonded p-orbitals that are 
only weakly hybridized, as in Fig. [TJ These networks 
exhibit a relatively small number of intrinsic over- and 
under-coordinated vertices. The latter are, in fact, re- 
sponsible for both the midgap states and the transition 
state configurations intrinsic to molecular transport in 
the quenched melts. In testing this hypothesis, we face 
the deeper question of the actual stability of aperiodic 
ppcr-bonded networks. Experiment shows (see below), 
that the enthalpy excess of a glass relative to the cor- 
responding crystal is less than the typical vibrational en- 
ergy, i.e., a small fraction - one percent or so - of the total 
bonding energy! In other words, despite their aperiod- 
icity, glasses are nearly defect-free structures, consistent 
with their bulk stability, both mechanical and thermo- 
dynamic. This observation is in conflict with a common 
view of glasses as a reconstructed - but otherwise arbi- 
trary - array of malcoordinated configurations and other 
local defects, such as vacancies. This common view would 
imply excess enthalpies of the order eV per several atoms, 
while avoiding to address the mechanism of transport in 
the melt. 




FIG. 1. Structure of rhombohedral arsenic as an example 
of a ppcr-bonded network. The solid lines denote regular co- 
valent bonds that connect the central atom with the near- 
est neighbors, bond length d. The dashed line denote the 
weaker, "secondary" bonds connecting the central atom with 
its next nearest neighbors, bond length d! . Angles /3 < 180° 
and a ^ 90° reflect the deviation from ideal octahedral co- 
ordination. These particular crystal fragment and view are 
based on Fig. 1 of Shang et al.— 



In the present and companion article,— we test the 
hereby proposed microscopic picture in two relatively 
separate stages. The present article is devoted to the 
first stage, in which we argue that ppcr-networks can rep- 
resent the quenched liquid and frozen glass forms of these 
substances in the first place. We will make a case that (a) 
such networks are stable against other types of bonding 
in a rather large class of amorphous compounds contain- 
ing elements from groups 15 and 16; and (b) despite their 
relative stability, aperiodic ppa networks are multiply de- 
generate, as are the actual liquids and glasses in question. 
The present work thus contains, to our knowledge, the 
first chemical bonding theory of a bulk glass. 

The article is organized as follows. In Section [TTJ we 
discuss in detail several key features of ppa networks, in- 
cluding their spatial non-uniformity and a hierarchy of 
bonding, from strictly covalent to weaker "secondary" to 
weaker yet van der Waals. Despite the presence of such 
weak bonds, the ppa networks are stable. To trace the 
origin of this stability, we formulate a structural model 
in Section inTl by which both periodic and aperiodic ppa- 
networks are symmetry-broken versions of a highly sym- 
metric, strongly bonded structure. This view is simi- 
lar but distinct from the common view of many elemen- 
tal solids as Peierls distorted simple-cubic lattices^ The 
symmetry breaking is driven by several competing inter- 
actions, including in particular sp-mixing; these, never- 
theless, are only strong enough to perturb but not qual- 
itatively modify the basic ppa bonding. We will verify 
that the resulting aperiodic pper-bonded lattice satisfies 
the three requirements for the existence of the topolog- 
ical states listed above. Importantly, this lattice will be 
argued to exhibit multiply degenerate configurations that 
differ by precise coordination of individual atoms, consis- 
tent with the possibility of activated transport in the cor- 
responding quenched melts. The precise degree of degen- 
eracy and the possibility of activated transport are both 
intimately related to the questions of the concentration of 
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the corresponding transition-state configurations in the 
melt and the accompanying electronic excitations. The 
latter questions are analyzed in the companion article.— 



II. ppa-BONDED SEMICONDUCTORS: THE 
ROLE OF THE SECONDARY ppa-INTERACTION 

The goal of this Section is to provide a detailed descrip- 
tion of ppa-networks in pnictides and chalcogenides, as 
arising from sigma-bonding between p-orbitals that are 
only weakly sp-hybridized. Such a description is neces- 
sitated by the lack of systematic comparative studies of 
the electronic properties of chalcogenide and pnictide al- 
loys, even though their structure itself has received much 
attentioni£i21 The ppcr-bonding emerges subject to com- 
petition from other types of local ordering. The presence 
of several competing types of local ordering in chalco- 
genides is evidenced by their broad range of structural 
and electronic properties, as could be seen by comparing 
e.g. As2Se3^S and GeSe2^ While the former exhibits a 
distorted octahedral coordination and well separated s 
and p bands, the latter displays tetrahedral ordering and 
overlapping sets of s and p orbitals. At the same time, 
the two substances exhibit opposite trends in terms of in- 
trinsic and light induced ESR responsei 21 ' 22 ' 30 ' 31 Such si- 
multaneous trends are particularly vivid in the Ge^Sei-^ 
series' 52 for 1/3 < x < 1/2, which exhibit coordination 
ranging from tetrahedral (smaller x) to distorted octa- 
hedral (larger x). In this series, the octahedral ordering 
seems to correlate with the separation between s and 
p band a 33 ' 34 and light induced ESR, and anti-correlate 
with the presence of unpaired spins 31 - and the glassform- 
ing ability.— Vice versa, the tetrahedral bonding exhibits 
the opposite trend. 

When sp-hybridization is weak, each atom exhibits 
a distorted octahedral coordination, as exemplified in 
Fig. [TJ Two or three nearest neighbors are situated at the 
distance of the regular covalent bond in an almost right- 
angled geometry. Opposite to each of these covalently 
bonded neighbors, there is an atom at a distance that 
exceeds the sum of the corresponding covalent radii, but 
is closer than the sum of the corresponding van der Waals 
radiii 3 ^ Crystalline As, Se, As2Se3, GeSe are typical ex- 
amples of this type of coordination. It is appropriate to 
think of crystals that exhibit distorted octahedral coor- 
dination not as fully covalently bonded, but as networks 
consisting partially of fully covalent bonds and weaker, 
closed-shell interactions. In the physics literature, it is 
customary to call the stronger bonds "front bonds" and 
to call the weaker bonds "back bonds."— ~— Because 
the back bonds formally correspond to closed-shell in- 
teractions, chemists call them "secondary" bonds^2r— 
or donor-acceptor interactions^ 4 - or, sometimes, hyper- 
valent or 3-center bonds, where the distinction is only 
quantitative, if any. 42 Importantly, the secondary bonds 
are stronger than van der Waals interaction and are direc- 
tional, similarly to their strictly-covalent counterparts^ 




FIG. 2. A fragment of the As2Se3 crystal. The secondary 
bond between Asl and Sel is shown with a thin solid line. The 
deviation from the strict octahedral coordination is reflected 
in different values of covalent and secondary bond strengths: 
t' a /t a — 0.44 < 1, a deviation from their strict alignment 
A/3 = 180° — P = 21°, and bond elongation relative to the 
sum of the covalent radii: Sd = +0.07 A, Sd' = +0.9 A. The 
notations are the same as in Fig. [1] 



We will see that, in fact, the distinction between the sec- 
ondary and covalent bonds in ppcr-networks is not sharp. 
A common example of the coexistence of covalent and 
secondary bonding is crystalline As2Se3, which consists 
of puckered layers of AsSe3 pyramids, see also Fig. 05] 
below, whereby the layers are only loosely bonded. The 
pyramids are made of the stronger, covalent bonds, while 
the secondary pper-interaction accounts for the rest of the 
intralayer bonding, see Fig. [5J The interlayer secondary 
bonds are even weaker; nevertheless they have been ar- 
gued to be as strong as 0.3 eV in As2Se3? 4 £ i.e. signifi- 
cantly stronger than a typical van der Waals-like bond. 

Despite thousands of documented instances of sec- 
ondary bonding and the similarity in its properties across 
a broad spectrum of compounds ) 41 ! 42 ' 44 both the mech- 
anism and quantitative description of this type of bond- 
ing still appear to be a subject of debate. 41 Without 
claiming full generality, here we will presume the ex- 
istence and universality of secondary bonds based on 
those myriad documented cases. The energetic and spa- 
tial characteristics of these bonds will be treated empir- 
ically, using tight-binding (TB) formalism. Such an ap- 
proach is controllable insofar as the substances in ques- 
tion are insulators or poor conductors since under these 
circumstances^ both localized, Wannier-like and the de- 
localized Mulliken-Hund orbitals form a complete set of 
electronic wave-functions. Consistent with the closed- 
shell character of the secondary bonding, the TB expres- 
sion for the corresponding binding energy contains ex- 
clusively interaction between occupied and unoccupied 
orbitals of the constituent molecules. Indeed, a one- 
electron Hamiltonian that describes the interaction of 
two molecules A and B can be written as a block ma- 
trix, where the ket corresponding to the wave function is 
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— tss tsp t a t' a ~tn ref. Structure 

Gel 1.70 2.36 2.56 nc , g ! i_ 0.67|diamond-Ge 

As 1.17 1.60 3.10 1.66 0.79 q-As 

Se 1.11 2.10 3.37 0.64 0.92 trigonal-Se 

TABLE I. Transfer integrals (in eV) in particular crys- 
talline forms of several elements often present in chalcogenide 
alloys.— All integrals are for the nearest neighbors, except 
t' a . In t ss and t sp , the subscripts indicate the constituent or- 
bitals. ta and t' a denote the transfer integrals for the covalent 
and secondary ppa bonds, and tn for the ppir interaction. A 
graphical summary of the transfer integral definitions is given 
in Fig.rrrja), see also p. 23 of Ref.^ 

a vertical stack of the kets pertaining to molecule A and 
B: 

Matrices T~L A , T-L B are the Hamiltonians of the isolated 
molecules A and B, while V contains the corresponding 
transfer integrals. According to the standard perturba- 
tive expression,— the binding energy of two closed-shell 
molecules A and B, E hind = E^ — E A ot — Ef ot reads: 

y n m n m ) n 171 

where ri A \^ A ) = E A \^ A ) and U B \^) = E*\^). La- 
bels "occ" ("unocc") denote summation over occupied 
(unoccupied) orbitals of the molecules A and B. 

Now, the precise geometry of the ppcr-network is sub- 
ject to several competing interactions, combined with the 
precise stoichiometry and other many body effects: the 
ppcr-interaction, sp-mixing, and pp7r-interaction, such as 
often found in conjugated polymers. (See Harrison^! for 
an introduction to tight-binding methods.) All these in- 
teractions have comparable strength as can be inferred 
from the values of the corresponding electron transfer 
integrals. Table Q] compiles the values of these transfer 
integrals for important representatives from groups four- 
teen, fifteen, and sixteen. Elements of these groups are 
of particular interest in the context of amorphous semi- 
conductors, because of comparable electronegativity and 
suitable valency, of course. The close magnitude of the 
listed competing interactions implies that the local order, 
which in turn is strongly affected by the stoichiometry, 
plays the crucial role in determining which interaction 
will ultimately dominate. 

One may list several complementary ways to establish 
the presence and significance of ppc-bonding. A ppa- 
network reveals itself structurally in a weak deviation 
from the octahedral coordination. In addition to a nearly 
right-angled geometry, the disparity between the lengths 
of the secondary and covalent bonds should be modest. 
In the latter case, the ratio of the corresponding trans- 
fer integrals, t'^/tu, is not too small, implying a relatively 



uniform, stable network. Furthermore, in view of a nearly 
universal relation —t^/tg — 1/4, Ref4£, a large enough 
value of t' a I t a automatically guarantees that the effect of 
ppir interactions on the geometry is small. On the other 
hand, when sp mixing is weak and little ss bonding is 
present, the top of the valence band consists primarily 
of p-orbitals, while the s and p subbands are relatively 
well separated. Indeed, consider for the sake of argument 
two identical centers, each having one s and p orbital 
and three electrons. The p-orbitals are aligned. Within 
the second order in the sp-mixing, the one-electron ener- 
gies of the four resulting molecular orbitals are given by 
e s =F t ss — t 2 sp [e p — e s ± (t a + t ss )] _1 for the ssa bond and 
e p ±t a + t 2 . p {e p ~e s ±(t cr + t ss )]~ 1 for the ppa bond. If the 
s and p orbitals are sufficiently separated in energy, it fol- 
lows automatically that (a) the centers are ppc-bonded 
and (b) the effect of the sp-mixing on the ppa trans- 
fer integral of the bond is small: t sp j \jt a {e v — e s ) < 1. 
A combination of the photoemission spectra with tight- 
binding calculations, using the known crystal structures, 
strongly suggests that exactly this type of bonding oc- 
curs in the crystals of several archetypal chalcogenides: 
As 2 S 3 , As 2 Se 3 , and As 2 Te 3 42 

In Table UH we compile data on the deviation from the 
ideal octahedral coordination and the corresponding ppa 
transfer integrals, in several distinct compositions and 
stoichiometries characteristic of common chalcogenide 
and pnictide alloys. One observes that a high value of 
t' a /t a is indeed characteristic of ppa-bonded materials, 
whereby the angular deviation from the ideal coordi- 
nation does not exceed 10°. Conversely, a large value 
of the t' a /t a ratio, alone, is often a good predictor of 
ppcr-networking. Note also a subtle, but nevertheless 
significant trend that in such a network, the stronger 
bonds are somewhat longer than the sum of the cova- 
lent radii. Furthermore, this deviation is the more signif- 
icant, the shorter - and hence stronger - are the secondary 
bonds. This anti-correlation is a telltale sign of "trans- 
influence"— in which the weaker bonded atom donates 
electrons into anti-bonding orbitals of the stronger bond, 
a common feature with secondary and donor-acceptor in- 
teractions. For the trans-influence to take place, it is 
essential that the counterpart covalent and secondary 
bonds be in a near linear geometry so that the anti- 
bonding orbital on the stronger bond overlap significantly 
with the bonding orbital on the secondary bond4^ Note 
that tight-binding descriptions are consistent with the 
trans-influence. Indeed, we show in Appendix B that the 
bond order— for the AB dimer from Eq. ^ is approxi- 
mately given by the expression: 

(A occ .B u „ occ A unocc B occ ^ . |2 
\ n m n m ) K^n 1 ~'m) 

The first (second) double sum is twice the occupation 
of the antibonding orbitals of molecule B (A), as do- 
nated by the bonding orbitals of A (B), consistent with 
the donor-acceptor nature of the secondary bond. Alter- 
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crystal and Ref. 
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t after 



strong ppa-secondary bonding/hypervalency 



As 2 Te 3 


50 


90.8 


0.1 


0.1-0.4 


0.4-1.2 


3-21 


.4-.9 


r-As 


51 


96.7 





0.12 


0.7 


16 


.65 


moderate ppa- 


secondary bonding 








o-As 


52 


97.0 





0.09 


0.9 


14 


.5 


GeSe 


53 


97.0 


0.7 


0.18 


0.9 


17 


.5 


t-Se at 5 GPa 


54 


104.7 





0.05 


0.8 


15 
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As2Se3 


55 


97.5 


0.3 


0.02-0.08 


0.9-1.4 


10-20 


.2-.4 


As 4 Se 4 


56 


98.4 


0.3 


-0.04+0.11 


1.1-1.5 


2-15 


.2-.4 


t-Se at GPa 


54 


103.1 





0.03 


1.1 


19 


.3 


r-Se 


57 


101.1 





0.02 


1.1 


9 


.3 


Br at 5 K 
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0.02 


1.0 


10 


.3 


Br at 250 K 
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0.01 


1.1 


10 


.3 


a-m-Se 


59 


105.6 





0.00 


Z 1-1 
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AS2O3 (As) 


00 


96.5 


1.1 


0.0-0.1 


1.2 


23 


.3 


AsBr 3 (As) 


61 


99.0 


1.1 


-0.1+0.1 


1.6 


11 


.2 



sp -bonding 



GeSe 2 (Ge) 


62 


109.4 


0.7 


-0.06-0.03 


> 1-3 


8 


< .3 


h-Ge 


63 


109.4 





-0.03 


1.6 





.2 


diamond- Ge 


64 


109.5 








2.2 


30 


.1 



TABLE II. Geometry (notations from Fig. [TJ and tight- 
binding parameters of representative compounds in crystalline 
form. Lower entries correspond to progressively weaker ppa- 
bonding. a is the average bond angle (hybridization), e is the 
half-difference in absolute electronegativities, Ad and Ad' are 
the deviations of covalent and secondary bond lengths from 
the sum of the covalent radii, A/3 = 180 — /3 is the devia- 
tion from a linear geometry, t'^/tc is the ratio of ppa integrals 
for covalent and secondary bonds. Crystallographic abbrevi- 
ations: h - hexagonal, m - monoclinic, o - orthorhombic, r 
- rhombohedral, t - trigonal. The computational details are 
given in Appendix lAl 



natively, according to Eq. (j4|), the bond order strongly 
correlates with the binding energy. Since the sum of the 
bond orders on a given atom equals to the atom's valency, 
stronger secondary bonds imply weaker counterpart co- 
valent bonds. 

Above said, we should point out that a large value 
of the t' a /t a ratio and the perfect octahedral coordina- 
tion, separately or together, do not guarantee that the 
ppcr-bonding is the main contributor to the lattice stabi- 
lization. An obvious counterexample is provided by ionic 
compounds with the rocksalt structure, which exhibit the 
perfect octahedral coordination. Incidentally, using sto- 
ichiometry to impose a (distorted) rocksalt structure is 
not guaranteed to produce a pper-network either. For in- 
stance, whereas GeSe does indeed exhibit distorted octa- 
hedral coordination, AsSe forms a molecular crystal com- 
posed of As4Se4 units, whose symmetry is incompatible 
with uniform octahedral coordination, despite relatively 
strong secondary bonding, see Table [Til and the Supple- 
mental Material*!^ Finally, other competing types of lo- 




FIG. 3. A dimer of Br2 molecules. The geometry is optimized 
at the MP2 level, using the program Firefly^ with aug-cc- 
pVTZ basis set and RHF wave-function. The gray shapes 
show the lowest energy molecular orbital consisting of the 
valence p orbitals of the Br2 molecules (61st MO, Table IVl in 
Appendix [Bl . see also the MO diagram in Fig. 1121 The two 
hues reflect the sign of the wave function. The computational 
details are provided in Appendix [Bl 



cal order are present in solids formed by the elements 
from groups 14-16. For instance, in GcSc2 the coordina- 
tion of Ge atoms is tetrahedral, see Table [IT] Elemen- 
tal phosphorus and sulfur at ambient conditions exist as 
molecular crystals made of tetraphosphorus P4 and octa- 
sulphur Sa , respectively, that show no signs of octahedral 
coordination. 

Apart from the stoichiometry and peculiar types of 
local ordering, the actual degree of stabilization of the 
pptr-network crucially depends on the strength of the sec- 
ondary bonds since they account for at least a half of the 
total bonds. We are not aware of systematic ab initio 
studies of the dependence of the strength of these bonds 
on the bond length and the deviation from the precise 
octahedral coordination, in a crystal. Such studies are 
understandably difficult, as the actual chalcogenide crys- 
tal structures exhibiting this type of bonding are very 
complicated. For instance, the unit cell of As2Se3 has 
20 atoms. 45 Despite these complications, it is possible to 
obtain an accurate estimate of the strength of ppa sec- 
ondary bonding in semiconductors by analyzing the sim- 
plest possible micro- and macro-molecular systems that 
exhibit this type of bonding, i.e. dimers of diatomic halo- 
gen molecules and halogen crystals respectively. Specif- 
ically, bromine is an appropriate example, because the 
elements of interest are located in periods 3 through 5. 
The ground state geometry of the bromine dimer is shown 
in Fig. [3] Here the ppa molecular orbital of the l.h.s 
molecule is mixed with the in-plane ppn molecular or- 
bitals of the r.h.s. molecule, as can be seen directly in 
the electronic density distribution in Fig. [3] Consistent 
with this mixing, the strength of the ppa secondary bond 
between the bromine molecules exceeds 0.3 eV, i.e., sig- 
nificantly more than expected of a typical van der Waals 
bond. According to the estimates on a variety of pseudo- 
dimer structures by Anderson et al^i, the strength of the 
secondary bonding, relative to the covalent bonding, de- 
creases toward the r.h.s. in each period. The above figure 
for the binding energy of the bromine dimer thus gives 
us a secure lower bound on the strength of a secondary 




lo'GPa pressure 



FIG. 4. Pressure dependence of the lengths of the covalent (d) 
and secondary (d!) bond in rhombohedral arsenic, after Silas 
et al.— At a critical pressure, of about 10 1 GPa, the lattice 
becomes simple cubic, while the two types of bonds become 
equivalent. In the inset, we replot the framed region with d 
as a function of d! and the pressure as a dummy parameter. 
The resulting dependence illustrates the trans-influence of the 
covalent and secondary bonds, c.f. Figs. 1 and 2 of Landrum 
and Hoffmann 42 

bond, consistent with the data in Table HU Furthermore, 
bromine crystals are comprised of layers in which Br2 
units are arranged in nearly the same geometry as in the 
ground state of an isolated pair of B12 molecules, see 
Supplementary Material*^ 



III. ppcr-NETWORKS AS SYMMETRY BROKEN 
STATES 

Here we propose a specific structural model of ppa- 
network formation, both periodic and aperiodic, and ar- 
gue that aperiodic pper-networks are consistent with (a) 
the structural degeneracy of the corresponding solid and 
(b) the restrictions on the magnitude of the spatial vari- 
ation of the electronic transfer integral from Eq. (TTJ) de- 
rived in Ref J4 

In liquids and glasses, the first coordination shell is de- 
termined by the stronger bonds and is very similar to the 
first coordination shell in the corresponding crystals^ - — 
Much less is known about the precise configurations of 
the weaker-bonded, next-nearest neighbors. A useful 
cue is provided by the observation that the crystalline 
photoemission spectra of several archetypal ppc-bonded 
chalcogenides - AS2S3, As2Se 3 , and As2Te ^ 49 ' 72 are very 
similar to their amorphous counterparts. In this and the 
following Sections, we argue that, indeed, the local in- 
teractions specific to pper-bonded glasses are of the same 
origin as in the corresponding crystals: In both cases, the 
local structures result from a symmetry-lowering transi- 
tion from the perfect octahedral coordination and thus 
are comparably stable. The main corollary of this infer- 
ence is that a supercooled liquid or quenched glass can 
be sufficiently stabilized by the pper-network alone. The 
key distinction between the crystal and glass is that ow- 
ing to its aperiodicity, the glass is necessarily structurally 
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FIG. 5. Parent structures of the crystals of (a) elemen- 
tal arsenic and (b) black phosphorus, after Burdett and 
McLarnan3 

degenerate. 

Let us begin with crystals. It is long 

appreciate d 36 ' 73 " — that the structures of many poly- 
morphs of elemental pnictogens, chalcogens, and 
halogens can be regarded as distorted simple-cubic (sc), 
with the exception of several of the lightest elements, 
such as nitrogen or oxygen, which form molecular 
crystals. Upon increasing pressure to several tens of 
GPa, the structures approach the ideal octahedral coor- 
dination, while phosphorus and arsenic actually exhibit 
a continuous transition to the simple cubic structure 
(Refsi^ 6 - and references therein). Specifically in arsenic, 
which is rhombohedral (A7) at ambient conditions, as 
the bond angles approach the right-angled geometry, the 
ratio of bond lengths of the nearest to the next nearest 
neighbor grows.^ 8 " These two bond types correspond to 
the covalent and secondary bonds respectively. In the 
vicinity of the transition, the covalent bond increases in 
length, while the secondary bond continues to shorten, 68 
see Fig. |4l Landrum and Hoffmann provide correlation 
data on thousands of pnictogen and chalcogen com- 
pounds that clearly demonstrate similar trans-influence 
between the covalent and secondary bonds 4^ see inset 
of Fig. |H The molecular fragments in question are of 
the type X-Q-X, where Q=Sb, Te and X = F, CI, Br, 
I. Note the correlation in the inset of Fig. @] pertains to 
valencies 3 (2) for Sb (Te). The combined view of Fig. [4] 
implies that the distinction between the covalent and 
secondary bonds is not sharp, but is subject to precise 
local coordination and/or bond tension. In the strict sc 
limit, both the secondary and covalent bonds become 
equivalent and should be regarded as fully covalent, 
albeit hypervalentj^ 3 - Similarly in halide crystals, each 
pprx-bonded layer transforms into a square lattice at 
a sufficiently high pressure (80 GPa for bromine 7 -^). 
Isovalent binary compounds A IV B VI transform into the 
simple cubic structure not only upon increasing pressure, 
but also temperature^ 

To rationalize these observations, we hereby pro- 
pose the following structural model, which draws heav- 
ily on Burdett and coworkers' view of the structure 
of rhombohedral arsenic, black phosphorus and other 
compounds i 78 ' 79 Place the atoms at the vertices of the 
cubic lattice, so that each atom is exactly octahcdrally 
coordinated and the p orbitals are aligned with the prin- 
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FIG. 6. Parent structure for (a) the GeTe crystal and (b) 
GeSe and GeS crystals, c.f. Fig. [S] Note these structures are 
by one horizontal layer taller than the repeat unit. 
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cipal axes. All atoms are linked, the links corresponding 
to bonds. Remove links so that each atom obeys the 
octet rule, while making sure the remaining bonds on 
each vertex are at 90 degrees, not 180. This procedure 
could be interpreted as adding electrons to a rocksalt- 
like compound while breaking bonds, whereby each filled 
antibonding orbital transforms into a lone pair of elec- 
trons pointing away from the remaining bonds^ As a 
result, each pnictogen and chalcogen, for instance, will be 
three- and two-coordinated respectively, whereby all links 
pointing from an atom are at right angles. We call the 
resulting lattice the "parent structure." Third, estimate 
the energy of the resulting parent structure, using a tight- 
binding Hamiltonian, while assuming that the transfer 
integrals are significant only for the linked atoms. Now, 
those bond-breaking patterns that have a particular low 
energy are special. One should expect that to these spe- 
cial structures, there correspond crystals of actual sub- 
stances that exhibit a distorted octahedral coordination, 
in which the covalent bonds will precisely correspond to 
the links, while the missing links correspond to secondary 
bonds or weaker, van der Waals interactions. For in- 
stance, Burdett and McLarnan have shown there are 36 
inequivalent ways to arrange three-coordinated atoms on 
the cubic lattice with a repeat unit of size 2x2x2. Two 
of the structures correspond to the lattices of black phos- 
phorus and rhombohedral arsenic^ In actual materials, 
both of these lattices consist of double layers that are 
buckled and mutually shifted, compared with the parent 
simple cubic structure. Other specific examples can be 
found in Refs j22t— It is understood that although the 
simple cubic lattice is a convenient parent structure for 
many compounds, it is by no means unique in this re- 
gard. For instance, Albright et al^l mention two addi- 
tional formal ways to obtain the arsenic structure, i.e. 
by adding two electrons per atom to wurzite ZnS or by 
puckering graphite sheets. Yet what distinguishes the sc- 
like parent structure is that, like the actual material, it 
is pper-bonded, whereas the orbitals in the wurzite and 
graphite structures are sp 3 and sp 2 hybridized respec- 
tively. We point out that IV- VI compounds that are iso- 
electronic with arsenic can be obtained from the parent 



FIG. 7. A parent structure for a Pn2Cri3 crystal, such as 
crystalline As2Se3 and AS2S3, c.f. Fig. E£b). 



structures of arsenic or phosphorus, see Fig. [6] Finally 
note that the above rules for bond placement, i.e. three- 
coordinated pnictogens and two-coordinated chalcogens 
with right angles between bonds, can be formally re- 
garded as a subcase of the 64-vertex model^ which is 
the 3D generalization of the venerable 6-vertex model 
of ice and 8-vertex model of anti-ferroelectrics* 8 ^ In the 
present case, 8 configurations on pnictogen vertices and 
12 configurations on chalcogen vertices have finite en- 
ergies, while the rest are infinitely costly. This analogy 
implies the proposed model is generalizable to more com- 
plicated coordinations by assigning finite energies to the 
latter. 

The present structural model can be formulated for 
stoichiometries that can not be arranged on the simple 
cubic lattice, except if one allows for vacancies. A specific 
example of particular relevance for this work is archetypal 
pnictogen-chalcogenides of stoichiometry Pn2Cli3, such 
as As2Se3, that can form both a glass and a crystal. (Pn 
= "pnictogen," Ch = "chalcogen.") In Fig. [7] we show 
that coordination-wise and symmetry-wise, the struc- 
ture of this compound consists of double layers similar 
to those in the black-phosphorus parent structure from 
Fig. Hfb). Note that by placing the vacancies in a partic- 
ular fashion, we achieve two things simultaneously: On 
the one hand, the As2Se3 stoichiometry is obeyed, and 
on the other hand, the octet rule on both pnictogens and 
chalcogens is satisfied. In drawing individual double lay- 
ers, we have used the structural model of Vanderbilt and 
Joannopoulos^ see also Ref4£ Note that based on the 
actual density of As2Se3, the bond length in the parent 
structure in Fig. [7] would have to be 2.8 A, a reasonable 
number consistent with the presence of the ppa network 
in the deformed structure. The presence of vacancies 
in parent structures should not be too surprising: The 
archetypal phase-change material Ge2Sb2Tes is known 
to exhibit a (metastable) distorted cubic structure with 
vacancies. 86,87 
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FIG. 8. The top view of a portion of the Pri2Ch3 parent 
structure from Fig. [7] The labels of the top layer are upside- 
down, to help distinguish it from the bottom double-layer and 
indicate that adjacent double layers are related by inversion. 
The thin dashed lines and adjacent numbers indicate the cor- 
responding secondary bonds and their length in the deformed 
structure, see Fig. [9] The tilted rectangle drawn with dashed 
lines indicates the unit cell. 



The parent structure in Fig. [7] is clearly not unique in 
that we could have placed the double layers in several 
distinct positions relative to each other. The specific, 
"homopolar" arrangement in Fig. [7J was chosen because 
the sets of close neighbors in this arrangement and in the 
actual deformed structure seem to exhibit the greatest 
overlap, see Figs. [5] and UU Nevertheless, several other 
mutual positions are possible, which exhibit comparably 
overlapping sets of close neighbors with the deformed 
structure, and, in addition, minimize the distance be- 
tween the vacancies in the parent structure better than 
the specific realization in Figs. [7J and [5] For instance, con- 
sider shifting the top layer in Fig. [5] "north" by one lat- 
tice spacing. Incidentally, one notices that the vacancies 
in the parent structure become "smeared" in the inter- 
layer space of the distorted structure, see Supplementary 
Material*^ We point out that it would be impossible to 
"merge" vacancies in the Pn2Cli3 stoichiometry between 
each two double-layers, if we attempted to use an arsenic- 
like structure as the parent structure from Fig. [SJa) in- 
stead of the black-phosphorus structure from Fig. EJb) . 
This notion is consistent with the stability of the latter 
structure in the actual material. At any rate, during the 
distortion, the distance between linked atoms decreases, 
resulting in strong, covalent bonds. Conversely, the sec- 
ondary bonds will result partially from the cleaved links 




FIG. 9. A side view of the actual As2Se3 structure. The wavy 
lines indicate inter-layer nearest neighbors in the parent struc- 
ture, except the A-3 link, which is through a vacancy. Note 
that the bonded intra-layer atoms are automatically nearest 
neighbors in the parent structure. The lengths of the links are 
only partially indicative of the actual bond length because the 
bonds are not parallel to the projection plane. 



(usually intralayer) or new contacts that formed in the 
deformed structure. Figs. [7J|T|] indicate that because the 
parent structure is not unique, there is no one-to-one cor- 
respondence between missing links in the parent struc- 
ture and the secondary bonds in the distorted structure. 

It is obvious that in the process of cleaving the bonds 
in the sc structure, symmetry was lowered resulting in a 
twice bigger unit cell along the pertinent directions. For 
example, in drawing the structure in Fig. [7J we could 
have split the original sc lattice in double-layers in six 
equivalent ways, i.e., two along each coordinate axis. 
Likewise, there are two equivalent ways to buckle each 
double layer, in each of the (1,1) and (1,-1) directions. 
Further symmetry lowering occurs when we place pnic- 
togens, chalcogens and vacancies at the lattice vertices. 
While the presence of symmetry breaking itself is hereby 
obvious, its mechanism appears to be subtler. Several 
workers have argued the symmetry-breaking transition 
that results in the structures of arsenic and phospho- 
rus is Peierls-like) 25 ' 69 ' 77 ' 88 which is a cooperative analog 
of the Jahn-Teller (JT) distortion^ i n ID or quasi- ID 
solids. The precise degree of Fermi surface nesting^ 
requisite for such a structural instability, seems however 
to be subject to a specific approximation employed . 25 ' 91 
The special significance of near- octahedral coordination 
and the resulting trans-influence between the correspond- 
ing covalent and secondary bonds, with regard to the 
Peierls metal-insulator transition, can be viewed from 
yet another angle: Alcock^S points out that compounds 
in which the two bond types display similar length show 
noticeable metallic luster. 

We thus conclude that one should generally regard 
the symmetry breaking of periodic sc parent structures 
as a second-order (or pseudo) Jahn-Teller distortion, 89 
whereby strict electronic degeneracy is not required. We 
further note that the solid-state analogs of the second- 
order Jahn-Teller (JT) effect are also well known, such 
as the dimerization transition in a /leieropolymer^ 2 - or in 
coupled homopolar chains, such as polyacenef^S Here, the 
polymer chain is unstable toward dimerization so long as 
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FIG. 10. The Peierls distorted (AsH 2 )„ chain. The links 
and gaps correspond to covalent bond (bond order 0.9) and 
secondary bonds (bond order 0.1) respectively. 



the gap is smaller than the coupling to the symmetry 
breaking perturbation, see Eq. ©. Electronic interac- 
tions, too, can lead to an effective one-particle gapJ^ Al- 
though second-order JT symmetry breaking, whenever 
present, is partially hampered by the lack of degeneracy, 
it is still driven by the very same mechanism as during a 
strict Jahn-Teller-Peierls distortion. 

We note that the question of the mechanism of the 
symmetry breaking is generally distinct from that of the 
interaction that determines the relaxed structure upon 
the symmetry breaking. According to Section [Til the lat- 
ter interaction in distorted-octahedral coordinated com- 
pounds is dominated by sp-mixing. Indeed, Seo and Hoff- 
mann point out^i that the distortion of the relaxed struc- 
tures away from the simple cubic structure is stronger for 
lighter elements, consistent with stronger sp 3 hybridiza- 
tion in those elements. For the perturbation caused by 
sp-mixing to also contribute to the symmetry breaking 
itself, it is essential that this perturbation can be made 
periodic with the inverse period commensurate with the 
electron-filling fraction in the undistortcd structure, as is 
the case (1/2 = 1/2) for trans-polyacetylene and arsenic. 
Otherwise, the Peierls-driven destabilization is severely 
weakened (Refi^, Chapter 2.6). 

Specifically i/iree-coordinated lattices with right angles 
between bonds have a very special property that make 
them additionally unstable toward symmetry lowering, 
even if these lattices are aperiodic. Any such lattice 78 
can be thought of as made of linear chains, each of which 
consists of white and black segments of equal length in 
strict alternation. The junctions between adjacent seg- 
ments correspond to the lattice vertices, while the white 
and black segments themselves correspond to no-link and 
link respectively. In the case of strict octahedral coordi- 
nation, i.e., no mixing between distinct chains, each such 
chain can be thought of as a result of a Peierls distor- 
tion of a chain of equidistant atoms. As a result, even 
aperiodic parent structures are expected to be relatively 
(meta-)stable, not only the strictly periodic structures of 
arsenic or black phosphorus. We illustrate the symme- 
try breaking in an individual linear ppa system in the 
presence of significant sp-mixing, which will be also of 
use later. Shown in Fig. [TU] is a dimerized linear chain 
of hydrogen-passivated arsenics, (AsH2)„. The details of 
the electronic structure and geometry determination for 
the chain are provided in Appendix [C] Despite the pres- 
ence of sp-mixing and other interactions, the chain ex- 
hibits a clear pper-character: (see Fig. Q]for the notations) 



The As- As bond length are d = 2.48 A and d' = 3.0 A for 
the covalent and secondary bonds respectively. The bond 
angles are anAsH = 97°, /?a s AsAs = 150°. The secondary 
bonding is significant as witnessed by the value of the 
corresponding transfer integral: t' ~ 2.3 eV, i.e. almost 
a half of that for the covalent counterpart: t ~ 4.9 eV. 
The competing interactions, i.e., ss, sp-mixing, and ppix 
are significantly smaller (see Appendix [Cjl . implying the 
bonding is indeed ppa. The sp-transfer integral is as 
large as 57% of the ppa transfer integral, consistent with 
the earlier statement that the sp-mixing is the next lead- 
ing contributor to the geometry of the symmetry broken 
state, after the ppa interaction itself. As expected of a 
Peierls insulator, the chain exhibits a perfect bond alter- 
nation pattern. Lastly note that even though individual 
chains that make up aperiodic 3D parent structures un- 
dergo Peierls transitions, the symmetry breaking for the 
3D structures themselves is generally not quasi-one di- 
mensional and thus is not Peierls, in contrast with peri- 
odic systems, such as considered by Burdett and others. 

The parent structures with three links per atom plus 
vacancies, if any, see Figs. [5JT1 can serve as basic mod- 
els for compounds with distorted octahedral coordi- 
nation that consist of two- and three- valent elements. 
These structures also apply to compounds containing 
elements from groups 14, such as GeTe, GeSe from 
Fig. E] or archetypal phase-change alloys Ge2Sb2Tes and 
GeSb 2 Te4i^ Hereby each A IV B VI pair is isolectronic with 
a pair of pnictogens, implying these atoms are three- 
coordinated. The rest of the atoms are in the Pn2Cli3 
stoichiometry, and so the rules leading to the parent 
structure in Fig. [7] apply to these atoms. Even in the 
absence of A IV B VI pairing, symmetry breaking schemes 
can be proposed for substances where the number of four- 
coordinated atoms is large, as in GeP and TII^S or j3- 
tin,— and similarly for five-coordinated atoms.— Oth- 
erwise, an atom with coordination 4 or higher can be 
considered as a defect in a lattice of three-coordinated 
vertices. We will return to this important point in the 
companion articled 

The following picture of structure-formation in ppa- 
networked materials thus emerges from the above con- 
siderations: These materials may be thought of as 
symmetry-broken versions of a simple cubic structure. 
The symmetry breaking is an interplay of several types of 
perturbation: (a) Peierls-instability proper for extended 
near-linear chains; (b) cooperative second-order Jahn- 
Teller distortion that results from mixing of the p-orbitals 
with other orbitals, mostly s, and from electronic inter- 
actions; (c) steric effects due to vacancies, if the latter 
must be present owing to stoichiometry, as in Fig. [7j and 
(d) other coordination variations, as in the case of el- 
ements from groups 14 and lower. Even though these 
perturbations are strong enough to break the symmetry, 
they are still perturbations, so that the resulting bond- 
ing is primarily ppa. This statement can be quantified 
by comparing the strengths of the corresponding transfer 
integrals in the deformed structure, see Table Q] and the 
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(AsH 2 )„ example above. 

Now, one must recognize that upon geometric opti- 
mization, the lattice will generally be aperiodic, even 
when the link-breaking pattern in the parent structure 
is itself strictly periodic, let alone if we arranged the dis- 
tinct species or vacancies aperiodically or with a period 
incommensurate with the period prescribed by the elec- 
tron filling fraction. Let us now examine how such ape- 
riodic lattices maintain the ppa character, and hence the 
stability with respect to other types of ordering, while, 
at the same time, allowing for molecular transport. 

For concreteness, let us consider a specific prescription 
to deform the parent sc structure. For each atom, con- 
sider the links in the parent structure as vectors that 
start on the atom. Move each atom by a small distance 
in the direction which is the sum of its own vectors. For 
instance, if an atom has three links: (1, 0, 0), (0, -1, 0), 
and (0, 0, -1), move it in the direction (1, -1, -1). Analo- 
gously for a two coordinated atom, the displacement will 
be in the plane containing the two links. Now turn on 
the interaction in the form of non-zero transfer integrals, 
such as listed in Tables Q] and [TTJ Let the lattice relax, 
subject of course to the Coulomb repulsion between the 
ionic cores. As already mentioned, there is generally no 
one-to-one correspondence between the parent and the 
distorted structure: this implies there are multiple re- 
laxed structures and hence multiple metastable minima 
on the total energy surface. In discussing Figs. [7]|9l wc 
have mentioned that distinct parent structures for the 
Pn 2 Ch 3 stoichiometry can be obtained by shifting the 
double layers relative to each other. Because such a shift 
incurs bond breaking, these alternative parent structures 
are separated by barriers. It is understood that since the 
parent structure is generally aperiodic^ Figs. EHH] apply 
only to a small fragment of such an aperiodic structure. 
Now, since the distinct parent structures are separated 
by an energy barrier, at least one of them should be sep- 
arated by a barrier from the actual deformed structure, 
implying a presence of additional, metastable minima. 
When such metastable minima are few, the global po- 
tential energy minimum is easily accessible, as is proba- 
bly the case for the periodic parent structure of arsenic, 
see Fig. EJa). Elemental arsenic is, in fact, a poor glass- 
former. (Other distinct three-coordinated parent struc- 
tures exist ^ but most of them are energetically costly 
To summarize, the existence of distinct parent structures 
with shifted atoms is crucial for the present structural 
model to be consistent with the presence of molecular 
transport in ppcr-networks. 

When aperiodic, the distorted lattice will exhibit two 
key features: First, if the substance forms a periodic crys- 
tal, it will be lower in energy than any aperiodic coun- 
terpart. Second, because aperiodic structures are not 
unique, the lattice will be highly degenerate, as just dis- 
cussed. We can also understand the emergence of the 
degeneracy thermodynamically: Suppose the substance 
can crystallize. A mechanically stable aperiodic struc- 
ture, on the one hand, has a much lower symmetry than 



the crystal. On the other hand, the aperiodic structure 
corresponds to a higher energy and hence higher temper- 
ature. The only way to reconcile these two conflicting 
notions is to recognize that there should be a thermo- 
dynamically large number of nearly degenerate aperiodic 
structures separated by finite barriers. By virtue of bar- 
rier crossing events, the lattice is able to restore the full 
translational symmetry at long enough times; the latter 
symmetry is higher than that in a mechanically stable 
crystal. The lattice therefore corresponds to a liquid in 
the activated transport regime, if steady stated or to an 
aging glass, if below the glass transition^ 

The view of quenched liquids and frozen glasses as bro- 
ken symmetry phases is supported by an independent 
argument: According to Fig. [4j the symmetry broken 
phase corresponds to a lower pressure. Bevzenko and 
Lubchenko^ have shown that a covalently-bonded equi- 
librium melt can be regarded as a high symmetry lat- 
tice that has been sufficiently dilated and then allowed 
to relax into one of the many available aperiodic config- 
urations. Now, are the predicted structural degeneracy 
of emergent aperiodic ppc-networks and the barriers for 
activated reconfigurations consistent with the configura- 
tional entropy of actual materials? The contiguity be- 
tween covalent and secondary bonding, as illustrated in 
Fig. |H suggests that ppcr-networks support atomic mo- 
tions that do not involve bond-breaking but only a grad- 
ual change in the coordination, and as such, may be ther- 
mally accessible. These special atomic motions will be 
discussed in detail in the follow-up articled 

Available structural data are consistent with the preva- 
lence of ppu-bonding in vitreous chalcogenides with 
the stoichiometries in question. According to several 
studies / 1 ! 96 the nearest neighbor bond lengths in amor- 
phous arsenic chalcogenides are essentially identical on 
the average to those in their crystalline counterparts, 
but have a somewhat broader distribution. In view of 
the argued presence of trans-influence between the cova- 
lent and secondary bonds in these compounds, we con- 
clude that the secondary bonds in the vitreous samples 
are of strength comparable to those in the corresponding 
crystals, again subject to a broader distribution. Conse- 
quently, based on the applicability of the same general 
mechanism of ppc-network formation and the compara- 
ble bonding strength, we conclude both periodic and ape- 
riodic lattices exhibit the same type of bonding. Still, for 
the sake of the argument, suppose on the contrary that 
the bonding is dominated not by the ppa interaction, but 
by its leading competitor, i.e., sp-mixing, thus resulting 
in a predominantly tetrahedral coordination. The ra- 
tio of the filling fractions of the diamond and As2Se3 
lattices is approximately 1.18. At the same time, the 
densitie a 55 ' 96 ! 97 of the amorphous and crystalline com- 
pounds differ by less, i.e. is 4.57 vs. 4.81-5.01 g/cm 3 for 
As2Se3 and 3.19 vs. 3.46 g/cm 3 for AS2S3, implying that 
only a small fraction of atoms in these glasses, if any, 
might be regarded as tetrahedrally coordinated. Such 
a "defect" is analogous to a small region occupied by an 
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interface between two distinct lattices, which incurs a sig- 
nificant free energy cost. The companion article^ shows 
that a mechanism for a reversible formation of such de- 
fected configurations arises naturally in the present struc- 
tural model. 



IV. CONCLUDING REMARKS 

The main goal of this article was to establish the 
mechanism of bonding in semiconducting pnictogen- and 
chalcogen-containing quenched melts and frozen glasses. 
Representative substances are archetypal glassformers, 
such as As2Se3 and similar materials whose crystalline 
forms can be directly argued to exhibit ppc-bonding, 
based on their known structures, measured electronic 
density of states, and electronic structure calculations. 
We have formulated a structural model, by which both 
the crystalline and vitreous materials are seen to form by 
the same general mechanism, i.e., by symmetry-lowering 
and distortion of a high-symmetry structure defined on 
the simple cubic lattice. By combining this model with 
the limited structural data on the vitreous counterparts 
of those listed materials and similar compounds, we have 
argued the glasses are also ppa-bonded. 

Lowering of the symmetry by breaking the bonds in 
the fully connected simple cubic structure can be under- 
stood as lowering of the lattice dimensionality, similarly 
to what is seen in the (AsEb)™ chain in Fig. [TU1 which 
is a linear array of relatively weakly bonded dimers, or 
to what one finds in the parent structures from Figs. [3} 
[7J which consist of double-layers, possibly accompanied 
by further symmetry lowering. Papoian and Hoffmann 43 
have outlined general principles for the interrelation of 
dimensionality and deformation of high symmetry pe- 
riodic structures. These authors argue the bonding in 
the parent structures from Fig. [5] can be regarded as a 
Peierls distortion of a hypervalently bonded lattice built 
from electron-rich units, whereby the dimensionality of 
the lattice is lowered. In contrast, to preserve the di- 
mensionality, the electron-rich units comprising the hy- 
pervalent structure would have to be oxidized instead. 
In the former case, the deformed structure retains its 
original character, while in the latter case, the geometry 
is expected to change. For instance, a cubic Sb° lat- 
tice exhibits a (distorted) octahedral coordination, while 
oxidation would hypothetically result in a tetrahedrally 
coordinated Sb + lattice. The semiconducting alloys in 
question do exhibit relatively low variation in electroneg- 
ativity and, hence, the distorted octahedral coordination. 

The parent structures for substances at the focus of the 
present study, i.e., pnictogen and chalcogen containing al- 
loys enjoy a very special property: Each vertex in the par- 
ent structure is f/iree-coordinated, while the bonds are at 
near 90° angles. Under these special circumstances, the 
whole lattice can be thought of as composed of infinite 
linear chains in which bond and bond-gaps strictly alter- 
nate. In view of the weak interaction with the environ- 



ment, each of these chains can be thought of as a quasi- 
1D Peierls insulator with renormalized interactions. (In 
the eventual deformed structure, the chains are deformed 
and, likely, of rather limited length^) This observation 
allows one to extend the trends established by Papoian 
and Hoffmann in periodic crystals to aperiodic lattices. 
Since the symmetry can be restored partially or fully 
by high pressure, the argued view of glasses as lowered- 
symmetry versions of high-symmetry structures is sup- 
ported by independent arguments on (negative) pressure- 
driven glass transition in covalently bonded materials. 95 
We have also pointed out that the view of this type of 
symmetry breaking as a second-order, cooperative Jahn- 
Teller distortion may be equally justified. In addition to 
the Peierls instability proper, within each chain, the sym- 
metry breaking is also driven by local interactions that 
compete with the ppcr-bonding proper, primarily by the 
sp-mixing. 

The presence of competing interactions in these glass- 
forming materials is consistent with their structural 
degeneracy^ It appears that the strength of such com- 
peting interactions should satisfy certain restrictions. 
Specifically, if the sp-mixing is too strong, it destroys 
the ppa-bonding. Yet, if the mixing is sufficiently weak, 
the coordination can be made close to perfect octahedral, 
while decreasing the stability of the glass relative to the 
crystal. The latter trend is exploited in making opti- 
cal drives, using Ge2Sb2Te5 or similar compounds i 86 i 87 
Again consistent with the structural degeneracy of the 
ppcr-networks is the noted analogy between the bond- 
placement rules in the proposed structural model and the 
vertex models known to exhibit rich phase behaviori 83 i" 

The argued similarity of the bonding mechanisms in 
ppa-bonded crystals and glasses explains the puzzling 
stability of this important class of glasses. As mentioned 
in the Introduction, the enthalpy excess AH of the su- 
percooled liquid, relative to the corresponding crystal, is 
easy to estimate. It is directly related to the configu- 
rational entropy S c of the fluid, i.e. AH = S C T, save 
a small ambiguity stemming from possible differences in 
the vibrational entropies. The liquid configurational en- 
tropy varies between and about 2k b per bead, be- 
tween the glass transition^ and melting temperatures. 
The energy 0.8fcsT fl per bead amounts to less than 0.05 
eV per atom. The model thus allows one to reconcile two 
seemingly conflicting characteristics of quenched melts 
and frozen glasses. On the one hand, these materials 
exhibit remarkable thermodynamic and mechanical sta- 
bility, only slightly inferior to the corresponding crystals. 
On the other hand, these materials are also multiply de- 
generate thus allowing for molecular transport. 

Finally, this paper buttresses and clarifies some of the 
technical aspects of the original programme: The ppa- 
networks automatically satisfy all the necessary require- 
ments for the presence of the intrinsic midgap electronic 
states, as listed in the Introduction. First, a necessary 
requirement for covalent ppa bonding is that the elec- 
tronegativity variation e is not too large, Eq. @, lest 
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(b) 















FIG. 11. (a) A graphical summary of the definitions of the 
transfer integrals from Table [I] (b) An example of geometry 
used to compute transfer integrals t and t', see text. 





reference structure 


fcOV , A 





a-quartz-Ge02 crystal 


0.52 


Ge 


crystal 


1.225 


As 


see text 


1.20 


Se 


isolated helb*^ and rings 


1.17 


Br 


diatomic molecule 


1.140 


Te 


isolated helix and rings^i 


1.37 



TABLE III. Covalent radii. The three columns contain the 
species name, the reference structure, and the resulting ra- 
dius. 



the resulting bonding becomes predominantly ionic. Sec- 
ond, a high-symmetry ppa-bonded network is unstable 
toward Jahn- Teller distortion at each center, implying 
the deformed lattice shows an alternating pattern of bond 
saturation in the form of covalent and secondary bonds. 
Third, based on the stability of the pper-network with re- 
spect to perturbations in the form of competing chemical 
interaction, the secondary bonds are sufficiently strong, 
i.e. the t'/t ratio of the transfer integral of the secondary 
and covalent bonds is not too small, see Eq. ([TJ. 
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Appendix A: Detailed explanations of Tables [I] and 

he 

Tabic |U A graphical summary of the definitions of the 
transfer integrals is given in Fig. [TTT a'). see also p. 23 
of Ref. 37 By definition, the ppa integrals, i.e., t a and t' a 
are computed assuming the two p orbitals are aligned 
with the bond and hence depend only on the distance 
between the participating atoms. As a result, the transfer 
integrals t and t' from Eq. Q] are equal in value to t a and 
t' a respectively only when the three atoms in question 
are situated on a straight line. If 180 — (3 < 20°, as 
is the case for almost all compounds from Table Q] the 
difference between t and t a is less than a few percent, 
depending on the specific geometry and basis set. For 
instance, in the geometry from Fig. Illf b). we get t = 
(p\\ % \pl) = t a + (t 7T — t a ) cos 2 (/3/2), yielding a difference 
of 4% between t and t„ for /3 = 160° and t v /t a = -1/4, 
Ref4£ 

Table HU For several compounds, the structural data 
vary somewhat depending on the source, such as for 
As 4 Se 4 in Refs.-^^ and a-monoclinic Se, Refs^^i The 
resulting ambiguity should be kept in mind. 

2nd column, deviation from the right-angled geome- 
try. For AS2O3, GeSe2, and AsBr3 only As and Ge were 
considered as the central atoms. 



3rd column. The absolute electronegativity is the av- 
erage of the electron ionization and affinity energies, as 
found in Ref^21 

4th, 5th column. The covalent radii used are listed in 
Table IIIII These radii were determined using compounds 
exhibiting low variations in electronegativity, as perti- 
nent to the materials in question, except Ge02- For As, 
values found in the Iiteraturei22r-i2i differ within 0.02 A. 
We use the value from the middle of this range, which 
also happens to coincide with the result of interpolation 
across the sequence of Ge-As-Se-Br. 

6th column. Only the strongest secondary bonds are 
cited. To determine these for As4Se4 and Se, we have 
used diagrams that show the magnitudes of the atomic 
p orbitals of the nearby atoms on a sphere centered on a 
chosen atom with a radius equal to the covalent radius of 
that atom, see the Supplementary Material*^ For the 
tetrahedrally bonded materials, AS2O3, and AsBr3, the 
nearest neighbor in the direction opposite to the covalent 
bonded atoms is used. 

7th column. Accurate values of the tight-binding 
(TB) parameters are usually determined by obtaining the 
best fit to the electronic density of states for each spe- 
cific systems and are thus system-dependent. Our goal 
here is, instead, to highlight the generic trends of tight 
binding parameters that apply satisfactorily for bonds 
ranging from the strictly covalent to van der Waals, 
in as many distinct compounds as possible. Such a 
universal parametrization of one-electron transfer (res- 
onance) integrals is provided, for instance, by the PM6 
parametrizationiS 7 - of the MNDO methodji2£ which we 
use to estimate the t' cr /t (7 ratios for the compounds in 
Table |TTJ Despite the made approximations, we feel that 
the resulting potential ambiguity of the ratio of the ma- 
trix elements, t'^/ta, is not large. 

We can partially judge the reliability of the tight- 
binding parameters by comparing them to those obtained 
by other standard methods. In Table ITVl below, we com- 
pare Robertson's data 39 for transfer integrals in three 
elemental solids, as obtained by fitting the spectrum 
of the one-electron TB Hamiltonian (column 2) and by 
the chemical pscudopotential method (column 3) to the 
present TB parametrization (column 4). The Figure on 
the r.h.s. compares the result of the parametrization for 
elemental Ge with accurate calculations of Bernstein et 
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Ref. 39 


Ref. 39 


present 




TB fit 


pseud. 


work 


As 


0.54 


0.51 


0.62 


Se 


0.19 


0.26 


0.27 


Te 


0.33 


0.65 


0.74 




5 3.5 4.5 d„ 



TABLE IV. Comparison of the t'^/tc ratio, as resulting from 
the present tight-binding parametrization with compilations 
of Robertson 39 (table) and the calculations of Bernstein et 
31,109 (figure^ see text. 



MOs 


AOs 


energy range, eV 


-Bbinding, eV 


61,62,64,66,68,69 




(-14.90,-10.9) 


-0.284 


57, 58,59,60 


S 


(-30.3,-25.7) 


-0.195 


63, 65, 67, 70 


Vz 


(-14.0,-10.9) 


-0.075 


71, 72 (unocc.) 


& 5 Px ,y 


(-0.3, +0.2) 


+0.466 



TABLE V. Bonding contribution of the molecular orbitals 
(MOs) of the Br 2 -Br 2 dimer from Fig. [3] The quantity 
-Ebmding is the difference between the total MO energies of 
the dimer and those of the isolated molecules. 



a l ,109 Finally, the same TB parametrization was used for 
calculation on the dimer of Br2 molecules, see below. 



Br Br 



Br- 



-Br- 



-Br 



Br 



Appendix B: Dimer of bromine molecules 

In this appendix we use the electronic structure of Br2- 
Br2 dimer obtained by ab-initio calculations to estimate 
the contribution of the ppa interaction to the binding 
energy of the dimer and show that the resulting estimate 
is consistent with the conclusions of a simple molecular 
orbital theory. 

The geometry of the dimer (see Fig. [3]) is optimized by 
Firefly 66 program on MP2-level with aug-cc-pVTZ basis 
set and RHF wave- function. Although more accurate ap- 
proximations exist the present method has the advan- 
tage of simplicity while yielding the geometry consistent 
with that of the bromine crystal; it also yields results that 
compare well with accurate calculations for other halogen 
dimersAii For clarity, we have not corrected for the basis 
set superposition error, because here we are interested in 
the relative magnitudes of distinct contributions to the 
bonding, not the absolute value of the binding energy. 
The resulting ambiguity in the binding energy itself is 
not large anyway, about 15%, as can be estimated using 
the standard counterpoise method, Ref *H p. 714. 

The total binding energy of the dimer, 0.13 eV, consists 
of several contributions. For instance, the correlation 
(MP2) contribution is 0.2 eV. To extract the contribu- 
tion proper of the electrons occupying distinct molecular 
orbitals, we adhere to the following scheme: First, we 
subtract the uniform downshift of all MOs by 20 meV 
in the dimer, relative to the two isolated Br2 molecules. 
Next, we compile the contributions of all molecular or- 
bitals to the bonding, see Table [V] Here, the dimer is 
in the x-y plane. Because little sp-mixing is present, the 
MO's are naturally grouped into classes that consist pri- 
marily of either p or s orbitals. The p orbitals are further 
subdivided, according to their symmetry, into p XtV and p z 
orbitals. The former form the ppa bonds, while the lat- 
ter are out of the x-y plane and contribute little to the 
bonding, as is clear from Table [V] 

The magnitude of the binding energy due to the in- 
plane p-orbitals, 0.284 eV, can be rationalized using a 



£,ev 
o - 

-1 - 

-10 - 

-11 - 

-12 - 
-13 - 
-14 - 
-15 - 



-72 
-71 



-H- 



3* 



44-69 
—2-68 



4*- 



4i- 



-H- 



.62 
■61 



4i- 



FIG. 12. MO diagram of the Br 2 dimer from Fig. [3] that de- 
picts a subset of the MO's from Table [Vj along with pertinent 
MO's of the constituent Br 2 molecules. 



simple molecular orbital theory. Let us first consider only 
the eight valence p-orbitals lying in the plane of the dimer 
shown in Fig. [3] The ab initio diagram of the correspond- 
ing MO's is shown in Fig. [T^l Let us focus exclusively 
on the intramolecular ppa- and pp7r-interactions and the 
intermolecular ppc-interaction between two atoms con- 
nected by the dashed line in Fig. |3l In this approxima- 
tion, the p y orbitals do not interact with the p x orbitals, 
so that the secondary bonding is exclusively due to the 
p x orbitals. The specific realization of Hamiltonian Q 
for this bonding configuration reads 



u 



(e p t 

trr £ 



o o\ 

p t* o 

£<j €p ^7T 

\0 U e p J 



(Bl) 
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By Eq. ((4]), combined with 



^unocc — e p + t a , 



"occ t p - 1 - L n 



we obtain the dimcr binding energy: 

(OVV (B2) 
Generally, in the perturbative expression (j4]), the quan- 
tities V nm correspond to the transfer integrals of the sec- 
ondary bonds, t', while the denominators to the transfer 
integrals of the covalent bonds, t. The binding energy is 
thus second order in the t'/t ratio. 

Note that the result in Eq. (|B2[) also helps to partially 
assess the reliability of the TB parametrization from Ap- 
pendix [AJ Substituting the numerical values for those 
transfer integrals in Eq. (IB2|) yields the dimer bonding 
energy 0.3 eV, in good agreement with the ab-initio anal- 
ysis leading to Table |Vj A similar analysis can be used 
to show that the s-orbital bonding contribution from Ta- 
ble [V] results not from the ss interaction, but primarily 
from sp-mixing, also consistent with the small overlap of 
the s-orbitals. 

Finally, we outline derivation of the TB expression 
([5]) for the bond order of a closed-shell interaction. Us- 
ing Eq. (14.22) of Rcf.— , for the wave-functions of an 
AB dimer formed by a closed shell interaction, one can 
straightforwardly show that the interaction-induced cor- 
rection to the density matrix V for the dimer, in the 
second order, is given by 



A occ B unocc 

E E 



A uno cc B occ 

E E 



Vn 



l^>«l 



(B3) 



According to a standard definition of the bond order j££ 
(also used in MOPAC), the latter can be written in the 
chosen basis set as: 



A B 




(B4) 



leading to Eq. 

Note that the first (second) double sum in 
Eq. ([5]) is twice the occupation of the antibond- 
ing orbitals of the molecule B (A), as donated 
by the bonding orbitals of A (B). These occupa- 
tions are defined as 2^| / AB ^ occ Em" nocc I (^ m |Vv) | 2 and 

2E^ AB)o ° c Em n ° c ° I (V4|VV> I 2 respectively, where V* 
are MO's of the AB dimer. 



Appendix C: Arsenic chain 

In our model system, the dimerized (AsH2)„ chain 
fFig.riU]). the shorter and longer As- As bonds correspond 



to covalent and secondary bonds respectively; the As- 
H bonds correspond to covalent bonds in chalcogenide 
glasses. The calculations are performed on semiempirical 
level by MOPACi^ program with PM6 parametrization. 
On each arsenic, the axes are labelled as follows: The p x 
and p y orbitals are oriented toward the hydrogens (but 
not strictly along the x and y axis) , while p z orbitals are 
directed toward the neighboring arsenics, see the geome- 
try in Fig. ITTT b). Note these axes are defined locally on 
each individual arsenic; the optimized chain may or may 
not be linear on the average. 

We first check the semiempirical approximation 
against ab-initio calculations for n = 4, see Table IVII 
The ab initio results clearly indicate a distorted octahe- 
dral geometry and are consistent with secondary bond- 
ing between the terminal (AsH2)2 units: The ft (As-As- 
As) angle is very close to 180°, while the bond length, 
3.57 A is close but shorter than the sum of the reported 
van der Waals radii for As, i.e., between 3.7 and 3.9 A, 
depending on the conventions^ Now, judging from the 
angle /? between adjacent As- As bonds and H-As-H an- 
gle, the semiempirical method clearly overestimates the 
tendency toward sp 3 hybridization. At the same time, 
the semiempirical method underestimates the secondary 
bond length, suggesting the two errors in the resulting 
bond strength will partially cancel out. The discrepancy 
in the secondary bond lengths and the /? angles between 
the ab initio and PM6 methods should be considered 
a consequence of the specific parametrization adopted 
in MOPAC, which was optimized for covalent, not sec- 
ondary bonds, and many of which have a significant ionic 
component and/or coordinations distinct from distorted 
octahedral. Thus the TB-based analysis is not expected 
to be fully quantitative, unless ad hoc parametrized. 
Nevertheless, we have seen from the discussion of Eqs. ^ 
and (0 that TB analysis is qualitatively consistent with 
the closed-shell character of the secondary bonding and 
the trans-influence between the counterpart covalent and 
secondary bonds. Our results for geometric optimization 
of rhombohedral arsenic, using the PM6 parametrization, 
are consistent with these conclusions: Using a 96 atom 
supercell, we obtain d — 2.46 A (vs. experimental 2.52 
A), d! = 3.24 A (3.10 A), a = 99.3° (96.7°). Clearly, 
the PM6 parametrization is consistent with ppa bonding 
in rhombohedral As, and, furthermore, gives reasonable 
figures for the lattice parameters. 

According to the semi-empirical calculation, an infinite 
AsH2 chain is dimerized and has a zero net curvature, 
whereby d = 2.48 A, d! = 3.0A, a = 97°, (3 = 150°. De- 
spite the aforementioned bias toward sp-mixing, analysis 
of localized MOs in the semiempirical calculations shows 
that the bonding contribution of arsenic's s-orbitals is 
negligible: 0.1 (vs. 2 for a regular bond). Significant 
secondary bonding and the trans-influence between the 
covalent secondary bonds are clearly seen in the values 
of the bond order, which are equal to 0.9 and 0.1 for the 
covalent and secondary As- As bond respectively. 

To independently verify the dominant character of the 
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dAsAs ^AsAs P dAsH QHAsH QAsAsH 


E,eV 


PM6 


2.463 3.06 148 1.523 95.7 97.4 
1.523 95.6 96.9 


-0.17 


MP2 


2.483 3.57 177 1.525 91.7 92.6 
1.527 91.6 92.1 


-0.13 


Ref^ 


(2.441) 3.53 (180) (1.506) (93.9) (93.9) 


-0.10 



of the effective Hamiltonian: 



TABLE VI. The geometry of the dimer of AS2H4 molecules: 
semiempirical PM6 calculations, ab-initio RHF MP2 calcula- 
tions with acc-pVDZ basis set, and calculations from Ref.-^ 
The notations are according to Fig. [T] two entries per cell cor- 
respond to the inner (1st) and the outer (2nd) AsHb units, the 
energy is the dimer binding energy, the values in parentheses 
are not optimized. 



4>\1> 


Pz 


s 


Px 




renormalization 


sp 


Pz 


-4.9 


-0.6 


0.2 


0.1 


e-e = +1.2eV 


60% 


S 




-13. 


0.9 


-4.0 






Px 






-4.4 


-7.6 






SH 








-5.0 






Pf 


4.9(f) 


2.8 


-0.8 


0.6 


i-t= +0.7eV 


50% 






-0.7 


0.4 


-0.1 












1.3 


0.4 






pt 


2.3(t') 


-1.0 


0.7 


-0.2 


i - t' = -0.3eV 


35% 






0.0 


0.2 


-0.1 






p% b 






0.5 


0.2 







TABLE VII. Columns 2-5: Elements {<f>\H\ip} of the Fock 
matrix of an infinite (AsH2) n chain. The transfer integrals t 
and t' are indicated in brackets next to their numerical values. 
Column 6: The renormalization of the parameters of the one- 
particle ppcr-only Hamiltonian, as computed using Eq. (|C2[l . 
Column 7: The contribution to this renormalization of the s- 
orbitals. Labels "cb" and "sb" mean covalent and secondary 
bonds. All AO's belong to the arsenics except the hydrogen's 



ppcr-bonding, we demonstrate that the competing inter- 
actions are perturbations that mostly amount to a renor- 
malization of the pper-interaction. Indeed, the full one- 
particle Hamiltonian for a system plus the environment 
is a block matrix 



H = 



H 



sys 

v n 



(Cl) 



where the matrix % sys contains exclusively the on-site 
energies and transfer integrals for the system, U cnv for 
the environment. The (generally non-square) matrix 
V and its hermitian conjugate contain the system- 
environment transfer integrals. It is straightforward 
to show that, given an energy eigenvalue E for the 
full Hamiltonian in Eq. (IC1[) . the portion of the wave- 
function corresponding to the system is an eigenfunction 



ft — ft^y S 



v* (E-n e 



V 



(C2) 



with the same eigenvalue E. Indeed, substituting \tp) = 
((■0 sys |, (ifj cnv \V into H\ip) = E\ip), one gets 'H sys |V'sy S ) + 
V*|Vw) = ^l^sys) and (^env - E)\tp cm ) = -V|^ sys ), 
from which Eq. (|C2j) follows. 

Clearly, in the one-particle picture, the effect of the 
environment can be presented as an (energy-dependent) 
renormalization of the bare Hamiltonian of the system 
proper. One can use this systematic procedure to esti- 
mate the strength of both the intra-chain and environ- 
ment's perturbation to the ppa bonding within the chain. 
In doing so, below, we fix the energy E at the gap center. 

The contribution of s- or d-orbitals to the renormal- 
ization of the ppa transfer integrals can be estimated by 
using a perturbative expansion of the exact Eq. (IC2|) : 



Hi 



Us 



^ E 



E 

a/3 



(E - e a ){E - e f3 ) 



where the Latin indices label intra-chain p z orbitals and 
the Greek indexes label the rest of the orbitals. The small 
parameter t/(E — e) is indeed small, and the more so the 
further the orbital energy e from the gap center. 

The essential elements of the Fock matrix are listed in 
Table IVnl Note that d-orbitals are also included in PM6 
parametrization. Despite the large value of the perti- 
nent transfer integrals, the direct contribution of these 
orbitals to the renormalization is at most a few per- 
cent because the orbitals are almost empty. One infers 
from Table I VIII that the dominant contribution to As- 
As bonding stems from ppu-integrals. The pp7r-integrals 
are four times smaller than ppa, while the ss-integrals 
are essentially negligible. Although s-orbitals do not by 
themselves contribute significantly to the bonding in the 
chain, they provide the main contribution to the renor- 
malization of the ppa transfer integrals, specifically by 
weakening the secondary bonding. This contribution, 
provided in the last column in Table IVII1 turns out to 
be well approximated by the expressions: 



Ac 



At « 2t sp ^- atom 



A. 



(C4) 



within the error of 20% or less compared with the more 
accurate Eq. (|C3|) . Here 



(C5) 



where g sp is the Coulomb integral for s and p orbitals of 
the same atom (6 eV for As) and P sp is the correspond- 
ing entry of the density matrix, which is the measure 
of actual hybridization. If P sp > - as is the case for 
AsH2 chain - then the sp-mixing weakens the secondary 
bonding, the magnitude of the effect proportional to the 
hybridization strength. 
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FIG. 13. Voids in As2Se3 crystal. The graph is generated as follows: First, we compute the distance r m i n to the nearest atom 
in each point (£, 77, Q in the lattice coordinates. (The lattice is monoclinic.) Call the resulting function r m i n (£, r/, Q. Next, we 
make a projection of this function onto the plane of Fig. 9 of the main text, which is perpendicular to the £ axis, according to 
the following recipe. For each point (£,??), we vary £ to find the largest value of r m i n (£, Tj, £), call it r(£, rf), and assign a color 
to the point (£,??) so that redder hues correspond to larger values of r(£,rf), more violet hues to smaller values of r(£, rf). The 
function r(£,rj) varies in the [1.64, 2.52] A range, while the sum of the covalent radii is 2.4 A. As a result, intensely red regions 
mark voids in As2Se3 crystal. Atoms are marked by white balls. The largest void is distorted octahedral; it is centered at point 
(1/2, 1/2, 1/2) with Se atoms as vertices. 




FIG. 14. Fragment of a 



slice of bromine crystal. Compare the bond angles with those in 



Fig. 3 of the main text. 



20 




FIG. 15. Illustration ol distorted octahedral coordination of the arsenic atom in the position Asl of As4Se4 crystal. Draw a 
sphere around the Asl atom at its covalent radius (1.2 A). For each point, compute the logarithm of the sum of the radial 
parts of the p orbitals on the atoms within 4.5 A. Assign a color to the value of this function so that redder hues correspond 
to larger values and more violet hues to smaller values. The resulting map is Merkator-projected on a rectangular area, where 
the top edge corresponds to the north pole. The resulting map is shown in panel (a). In perfect octahedral coordination, red 
areas would be centered at the black dots. (The line connecting the poles coincides with a C3 axis of the octahedron.) In panel 
(b), the nearest neighbors are excluded, to specifically highlight back-bonding. The pertinent neighbors are indicated, together 
with their distance to the Asl atom, in Angstroms. 




FIG. 16. Distorted tetrahedral coordination of the selenium atom in the position Sel of As4Se4 crystal. See Fig. 
explanation. Black and white dots indicate octahedral and tetrahedral coordination respectively. 
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FIG. 17. Distorted octahedral coordination of the atoms in trigonal-Se crystal. See Figs. 1151 and [16] for explanation. 



